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1 Introduction 



Electromagnetic properties of nucleons are an important source of information on the internal 
structure of baryons. The success of the nonrelativistic quark model for the description of static 
characteristics (masses, magnetic moments, etc.) and the results from deeply inelastic lepton 
scattering are a clear indication for the three-quark structure of nucleons. 

In view of the difficulties of describing the nucleon as a relativistic three quark system rigor- 
ously, on the basis of QCD, many methods and models have been developed which implement 
important aspects of QCD at least partially. We mention a few. Some approaches [1],[2] are 
formulated using classical hadronic degrees of freedom (color singlets) without explicit reference 
to quarks and gluons. This holds in particular for chiral perturbation theory which addresses 
low energy scattering of hadrons [3]. In the Skyrmion model the baryon is considered to be a 
solitonic excitation of classical meson fields which are described by an effective chiral Lagrangian, 
see [4], e.g. The electromagnetic radii and magnetic moments of nucleons as well as the behavior 
of electromagnetic form factors up to 1 GeV^ are described in this way. 

In the context of quark models the nonrelativistic approximation is problematic even for the 
constituent-quark picture at low energies since the effective masses and the intrinsic momenta 
have the same order of magnitude. An attempt to implement relativistic invariance for the 
description of the electromagnetic properties of the nucleon is the covariant constituent quark 
model by Chung and Coester [5] which uses light front dynamics for the constituent quarks. 
In this framework the available data for < < Q GeV^ has been well described with two 
adjustable parameters (constituent quark mass and confinement scale). 

For large momentum transfers, » m^, perturbative QCD predicts that the magnetic 
formfactor behaves as Gm ~ ^IQ^ [6]- The use of QCD sum rules allows to extrapolate the 
electromagnetic form factors to low and moderate by incorporating local quark-hadron duality 
[7]. Recall that the calculation of hadronic form factors by sum rules in their original version 
[8] is restricted to intermediate momentum transfers and does not work in the infrared region 
due to power corrections l/Q^". Therefore, in order to calculate magnetic moments of nucleons 
QCD sum rules in the external field approach have been introduced [9]. For a general method of 
calculating the nucleon magnetic form factors at small see [10]. 

For applications at high energies the diquark model [11] has been proposed. In this model 
the proton is built from quarks and diquarks and the diquark is treated as a quasi-elementary 
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particle. Fits to the data for > 3 GeV^ [12] with few parameters have been obtained. 

A different hne of attack uses both quark and hadronic degrees of freedom. Chiral potential 
models (see, for example, [13] and references therein) fall into this class. They use an effective 
confining potential for the quarks and a quark pion interaction which preserves chiral symmetry. 
Incorporating the lowest-order pionic correction, the magnetic moments of the nucleon octet 
have been calculated [13]. Alternative approaches [14]-[17] start from effective Lagrangians which 
describe the transition of hadrons to their constituent quarks, combined with some assumptions 
on the behavior of quarks at low energies and on the shape of the hadron-quark form factor. This 
approach is connected fairly directly to the difficult issues of hadronization and quark confinement 
in QCD. 

In the present paper we extend the line of thought developed in [18] and [19] which use local 
hadron-quark or local hadron-quark-diquark vertices together with a confined quark propagator 
(no poles). It has been shown for these models that the electromagnetic form factors of nucleons 
tend to be somewhat above the data for < 1 GeV^. This might be a hint that the assump- 
tion of a local coupling of baryons with their constituents is not quite adequate. Introducing 
a nonlocal hadron-quark vertex may be a good way of effectively introducing other degreees of 
freedom (mesonic cloud, soft gluons, etc.). The Nambu-Jona-Lasinio models, NJL, with separa- 
ble interactions are candidates for this purpose. Many successful applications of such models for 
low-energy pion physics exist (see, e.g., [14, 15] and references therein). 

In [20, 21] the Lagrangian formulation of the NJL model with separable interaction has been 
given both for mesons and (for the first time) for baryons. The gauging of nonlocal interactions 
has been done by using the time-ordering P-exponent. The pion weak decay constant, the two- 
photon decay width, as well as the form factor of the 7*7r° — > 7 transition, the pion charge form 
factor, and the strong ttNN form factor have been calculated and good agreement with the data 
has been achieved with three parameters, the range parameters characterizing the size of mesons 
Am and baryons A^, and the constituent quark mass niq. In this paper we shall use the approach 
developed in [20, 21] in order to calculate the electromagnetic form factors of nucleons. First, 
we slightly modify the gauging of the nonlocal quark-hadron vertex by using a path-independent 
definition for the derivative of the time-ordering P-exponent [22, 23]. Second, we derive Feynman 
rules for the diagrams which describe the electromagnetic nucleon form factors. For simplicity, we 
use a Gaussian shape for the nucleon-quark vertex. The permissible range for the two adjustable 
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parameters, the Gaussian range A of the separable interaction and the constituent quark mass m^, 
is obtained by fitting the experimental values of the magnetic moments and the electromagnetic 
radii. The momentum dependence of the electromagnetic form factors is then calculated and 
their shape is found to be close to the data, for the parameters obtained from fitting the static 
properties. 

2 Model 

We start with a brief review of our approach. We consider the hadron as being a bound state 
of relativistic constituent quarks with masses niq [20, 21]. The transition of hadrons into con- 
stituent quarks and vice versa is described by the corresponding interaction Lagrangian, and 
the momentum distribution of the constituents by an effective relativistic vertex function. Its 
shape is chosen to guarantee ultraviolet convergence. At the same time the vertex function is a 
phenomenological description of the long distance QCD interactions between quarks and gluons. 
The interaction Lagrangian of nucleons with quarks is written as 

2 I 



C%\x) = N{x) J dy, j dy, j dy, 5 - ^^±^±1^^ p ^(y. - y,] 

X {dNJlviyi, y2, ys) + gNJN{yuy2, ys)} + h.c. 
= N{x) J dir j di.Fiil + ii) 



+ (/w4(^-2ei,^ + ei- V36,x + ei + V36)| + h.c. (i 



with J"^ and Jjj being the vector and tensor currents, respectively. The currents are symmetric 
under permutation of any two quarks 

4(yi,y2,y3) = rri''q'''{yi)q''Hy2)r2TC^,q^^iy,)e^^^^^\ (2) 

4(yi,y2,y3) = ra^'^7'9"Hyi)9"^(y2)r2fCa,.g'^^(y3)e"^"^"^ (3) 

The notation is as follows: {a} stands for the color indices, a^^ = |[7^,7''], and C = ^'^^'^ is the 
charge conjugation matrix, respectively. 

The notation for the spin-flavour structure of the nucleon currents is the same as in [9, 24]. 
In particular, the proton quark currents in terms of isotopic u and d flelds have the form 
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(4) 



Jp = cr^ ^ a u Ccr^^u e . [b) 

Both forms have been used separately in [18, 19] for studying the electromagnetic and strong 
properties of baryons. It was shown that the tensor current (5) is more suitable for the description 
of the data. For this reason we will use the tensor current (5) in the approach developed in this 
paper. 

The vertex form factor Fi^S^I + S^^) characterizes the distribution of u and d quarks inside the 
nucleon. The choice of variables in the vertex function Fi^S^I + S^^) implies the use of the center 
of mass frame in which the coordinates of quarks are expressed as 

yi=x-2ii y2 = a; + ^1 - ys = a; + ^1 + 

^ y2 + ys - 2yi ys - y2 

= — I — ^' = ^wr 

with (^1 and (^2 being the Jacobi coordinates of the quarks. 

The Fourier-transform of the vertex form factor is defined as 

To introduce the electromagnetic field into the nonlocal nucleon-quark Lagrangian (1), we will 
use the time-ordering P-exponential [20, 21]. The interaction of quarks with the electromagnetic 
field is introduced by the standard minimal substitution. 

Thus, we have the gauge invariant interaction Lagrangian 

= gNN{x) J dyi J dy2 J dy^8 - ^L+^_+Ii^ p |^_L ^(y^. _ y.^|^ 

yi y2 
X TiT expl^ieQ j dz^A^y^'{yi)exp(^ieQ j dz^A^y^'iy^) (7) 

X X 

ys 

X CT2 T2T exp (^leQ J dz^A"^ q"' (y3)e"'"'"' + h.c. 

X 

where Fi and r2 are spin matrices, and Q=diag(2/3,-l/3,-l/3) is the charge matrix of quarks. 



This interaction Lagrangian generates nonlocal vertices which couple nucleons, photons and 
quarks. In general several diagrams contribute. For example, the electromagnetic form factors 
of nucleons are described by the standard triangle diagram (Fig. la) and, additionally, by the 
bubble or contact diagrams (Fig. lb and Fig.lc). 




(a) (b) (c) 

Fig.l Nucleon Form Factor. Triangle diagram (a), bubble diagrams (b) and (c). 

At this stage the physical matrix elements describing the electromagnetic processes depend 
on the path which connects the lower and upper limits in the integral of the P-exponential. In the 
papers [20, 21] we used a straight line in our calculations. Here, we shall use a path-independent 
formalism which is based on the path-independent definition [22] of the derivative of the line 
integral. This formalism has been used in [23] for gauging nonlocal Lagrangians. We recall that 
the derivative of the line integral 



/(x,y,P) = J dz,A^{z), 
y 

where P is the path taken from y to x, may be defined as [22] 

d 

lim dXfj^— — I[x^y^P)= lim [/(x -|- c/x, y, P') — /(x, y, P)]. 

Here P' is obtained from the path P by adding the extension dx to the x end. 
Using the definition Fq.(8) leads to 

d 



^x^^ 



(9) 



This means that the derivative of the /(x, y,p) does not depend on the path used in the definition. 
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The method for deriving Feynman rules in the nonlocal chiral quark model with the elec- 
tromagnetic P-exponential has been developed in [23]. The crucial point is to calculate the 
expression 

F{-dl)e-^^''I{x,y,p) 

using Eq.(9). 

It is readily seen that 

F(-ae-^-/(x,y,P) = e-^^^F{-Vl)I{x,y,P) 

where T>1 = [d^ — ipY . Thus we need to calculate 

F{-Vl)I{x, y, P) = ^(-r yVl-Iix, y, P). 

One finds that 

Vll{x, y, P) = {d,A{x)) - 2i{pA{x)) - p'l{x, y, P) = L{A) - p'l{x, y, P) 

where L{A) = - 2i{pA{x)). 

Then 

{Vlfl{x,y,P) = {Vl-p')L{A) + {-pyi{x,y,P) 
{Vlfl{x,y,P) = {{Vir-Vy+p')L{A) + {-p'fl{x,y,P) 

n-l 



iVlriix,y,P) = Y.(T^lr'-\-p'fL{A) + {-p'ri{x,y,P) 

k=0 

1 

= nj dt[Vlt - /(I - t)r-'LiA) + i-p'riix, y, P) 



Finally, this leads to 

1 

Fi-Vl)Iix, y,P) = J dtF'[-Vlt + p'il- t)]LiA) + Fip')Iix, y, P) 



1 X 

= f dqA,{q)\^(q^ + 2p^y-'^^ I dtF'[{p + q)H + p^l - t)] + F{p^) j dz^e-'^^^Y 

y 

We observe that the last term vanishes due to the delta function 8{x — y) in the Lagrangian. 
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For baryons, one has to calculate 

F{-dl - d^)e-'P'^'-'P'^' j dz^A^{z). 

y 

By analogy we get 

(-zai)(aig^ + 2pt)F, + {-za2){a2q^ + 2p^)F2 

y 

with Fi and F2 being defined as 

{P2 + a2qf + t{pi + (iiqf + (1 - t)pl 





1 

F2 = I dtF' 


With this preparation the Feynman rules for the ^NN vertex function may be derived. The 
details are given in Appendix A for the evaluation of the vertex in Fig. 2 which connects a 
nucleon, a photon and three quarks, while appendix B gives the connection to the nucleon mass 
operator of Fig. 3. 



dqAf,{q)e~ 



«>lCl-«>2C2 J „--iq(x+aiCi+a2C2) 



Fi 



dtF' 



Pl+t{p2 + a2qf + {l-t)p: 



(10) 




Fig. 2. Nucleon photon three quark vertex. Fig. 3. Nucleon mass operator. 

Now we turn to a discussion of the model aspects of our approach and the parameters involved. 
First, there is the vertex function in the Lagrangian (1). This is a free function except that it 
should make the Feynman diagrams ultraviolet finite, as we have mentioned before. In the papers 
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[20, 21] we have found that the basic physical observables of pion low-energy physics depend only 
weakly on the choice of the vertex functions considered. Therefore, we choose in this paper a 
Gaussian vertex function for simplicity. In Minkowski space we write 



where Ajv is the Gaussian range parameter which is related to the size of the nucleon. 

The second input is the quark propagator. In the present paper we shall use the standard 
form, corresponding to a free fermionic field with mass niq 



Of course, the deficiency of such a choice is lack of confinement. This could be corrected by 
changing the analytic properties of the propagator. We leave that to a future study. For the time 
being we shall avoid the appearance of unphysical imaginary parts in the Feynman diagrams by 
restricting the calculations to mjv < 3mg. 

Thus, there are two adjustable parameters in our model: the constituent quark mass niq and 
the range parameter Ajv- We shall determine them by fitting the electromagnetic properties of 
nucleons. 

3 Nucleon Mass Operator and the jNN Vertex Function 

Within the relativistic quark model described above we shall evaluate the magnetic moments, 
the charge radii and the behavior of the electric and magnetic form factors for < < 1 GeV^. 

It is convenient to separate the contributions from the triangle diagram (Fig. la) and the 
bubble diagrams (Fig. lb and Ic) as follows 



which are orthogonal to the photon momentum g^A^(p,p') = 0. Thus the Ward-Takahashi 
identity following from (12) reads 








(13) 
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In particular, if g = and both nucleons are on the mass shell, one finds 

d 

The expressions for the orthogonal parts of the vertex functions are 



Ka{p,p') = malY 



{2^)H J {2^)H 
F(12[A;2 + k'^ + kk']) F(12[A;2 + k'^ + kk'] + V - 12A;g) 



36gZ^ (A;, p, p')J\pp{k, k',p) + SGQ-f.Z"^ (A;, p, p'h,Ilss{k, k\ p) 



+ rQra"P^iat{k,P,p'y"'l,Wr{k,k',p) 



{2t^)H J (27 

F(12[A;2 + k'^ + kk'] + 4g2 - 12A;g) - F(12[A;2 + k'^ + kk']) 

— 3kq 

X 36^5', (^k + npp(A;, A;', p) + 86^756', (^A; + -f^Iissik, A;', p) 
+ TQTa"^-f,Sg[k + ^a^^n^^^'^^lA;, A;',p) 
+ (p p') 



where = (^fc + - f j 7;^^, (^fc + § j , Q =diag(2/3,-l/3), f^f =diag(0,l), 
7u = 7a* - 9a* K = - Qf^k ■ q/q^ (7-^ • g = and A;-^ • g = 0). 



npp(A;, A;',p) = tr ^fsSgi^k' + 0755'g(^A; + ^' ~ 



n55(A;,A;',p) = tr 



n:^^^'''"(A;,A;',p) = tr 



P 
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The nucleon-three-quark coupling gjf is calculated from the compositeness condition (see [16]) 
which means that the renormalization constant of the nucleon wave function is equal to zero, 
Zjv = 1 — S^(mjv) = 0, with Ejv being the nucleon mass operator 

E.(.) = 216(.^)7^/|^F^(l2[.^ + .'^ + (15) 



The calculation of the electromagnetic vertex functions and of the nucleon mass operators is 
given in Appendix B. 

4 Numerical Results and Discussion 

We start with a summary of well-known definitions for the electromagnetic properties of the 
nucleon. The vertex functions in terms of Dirac pj^''^'^ and Pauli Fp'^"^ form factors are 

M, = N\p')K^^{p,p')N^{p), (16) 

where N{p) and N{p') are the spinors corresponding to initial and final nucleon states with 
momentum p and p\ respectively, and 

N = {pn), ^ = ( n ) ^^^^ 

The electric and magnetic (Sachs) form factors are obtained from 

• Electric Eorm Eactors G^r\f) = F^j^'''\q^) + ^F^f'^Vg^) 

. Magnetic Eorm Eactors Cfe'^H?') = F^^'''\q^) + ^^^'"^(g^) 
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At zero recoil, = 0, these form factors satisfy the following normalization conditions 



am = 1, Glio) = 

where fip and /i„ are the magnetic moments of the proton and the neutron, respectively. 
The charge radii of the nucleons are given by 
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g2=0 



Fm + 



1 



4m% 



:is) 



dGUq' 



E- 



dq^ 



6 



g2=0 



F-(0) + 



4m%. 



(19) 



6 



dGUq') 



M- 



GU^) dq^ 



6 

/ijV 



F^'(O) + Fj^'iO) 



Nit 



(20) 



We shall work in the limit of isospin invariance where the masses of u and d quarks are 
equal, m„ = = niq. In a first step we set the parameters of our model niq and Ajv to the 
well measured static properties of the nucleons, to the magnetic moments and to the charge 
radii [25, 26]. In a second step we shall explore the constraints from the shape of the nucleon 
electromagnetic form factors for the range <1 GeV^. 

In Table 1 the results for the static properties are given for various sets of parameters niq and 
Ajv- As was mentioned above, the constituent mass niq is chosen to be more than one-third of 
the nucleon mass [niq > niN/3) to avoid an imaginary part in the Feynman diagrams. Its value 
is varied in the region from 315 MeV to 500 MeV. One of the important constraints is the ratio 
of magnetic moments of proton and neutron R = which is measured with high accuracy 

and equal to R = 1.46. Historically, the first success in describing this value has been achieved 
in a naive nonrelativistic SU{6) quark model. The value Rsu(6) = 1-5 [27] is already close to 
experiment. It is seen from Table 1 that in our approach this value varies in the range i?=1.50- 
1.57 when the quark mass is smaller than 400 MeV, and starts to be less than 1.5 for niq >400 
MeV. The electric neutron radius squared, < >^, is also sensitive to the choice of the quark 
mass. In this case as well the best agreement with experiment is achieved for niq >400 MeV. 
As can be seen from (19) the neutron charge radius is a combination of the neutron Dirac form 
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Table 1. Static Nucleon Properties 



niq, Ajv 
in GeV 


jjjrp 




\Ep.\ 


p 

(fm) 


< >'% 
(fm2) 


p 

(fm) 


' M 

(fm) 


0.315, 2.40 


2.78 


-1.77 


1.57 


0.91 


-0.171 


0.84 


0.85 


0.320, 2.25 


2.79 


-1.79 


1.56 


0.88 


-0.169 


0.80 


0.81 


0.325, 2.12 


2.81 


-1.81 


1.55 


0.88 


-0.166 


0.79 


0.80 


0.330, 2.00 


2.82 


-1.82 


1.55 


0.88 


-0.162 


0.79 


0.80 


0.335, 1.90 


2.82 


-1.84 


1.53 


0.88 


-0.159 


0.79 


0.80 


0.340, 1.80 


2.82 


-1.85 


1.52 


0.89 


-0.155 


0.80 


0.80 


0.350, 1.70 


2.82 


-1.85 


1.52 


0.89 


-0.152 


0.80 


0.80 


0.360, 1.60 


2.83 


-1.86 


1.52 


0.89 


-0.148 


0.81 


0.81 


0.370, 1.50 


2.82 


-1.87 


1.51 


0.90 


-0.144 


0.82 


0.82 


0.380, 1.40 


2.82 


-1.87 


1.51 


0.91 


-0.139 


0.83 


0.83 


0.390, 1.35 


2.81 


-1.87 


1.50 


0.92 


-0.137 


0.84 


0.84 


0.400, 1.30 


2.81 


-1.87 


1.50 


0.92 


-0.135 


0.84 


0.84 


0.410, 1.30 


2.79 


-1.86 


1.50 


0.91 


-0.134 


0.83 


0.83 


0.420, 1.25 


2.79 


-1.86 


1.50 


0.92 


-0.132 


0.84 


0.84 


0.430, 1.25 


2.77 


-1.86 


1.49 


0.91 


-0.132 


0.83 


0.83 


0.440, 1.25 


2.76 


-1.85 


1.49 


0.90 


-0.131 


0.82 


0.82 


0.450, 1.25 


2.76 


-1.85 


1.49 


0.89 


-0.131 


0.81 


0.81 


0.460, 1.20 


2.75 


-1.85 


1.49 


0.90 


-0.129 


0.82 


0.82 


0.470, 1.20 


2.74 


-1.84 


1.49 


0.89 


-0.129 


0.82 


0.81 


0.480, 1.15 


2.74 


-1.84 


1.49 


0.90 


-0.127 


0.83 


0.82 


0.490, 1.15 


2.73 


-1.84 


1.48 


0.90 


-0.127 


0.82 


0.82 


0.500, 1.15 


2.73 


-1.84 


1.48 


0.89 


-0.127 


0.82 


0.81 


Exp. 


2.79 


-1.91 


1.46 


0.86±0.01 


-0.119±0.004 


0.86±0.06 


0.88±0.07 
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factor F^'[0) and the ratio Sfin/'^'m-N- Using the experimental value for the neutron magnetic 
moment the last term almost reproduces the experimental value for < >^. This means that 
the contribution coming from F^'[0) should be negligibly small (Fg'(O) ~ 0.008). One can see 
from Table 2 that this value, disregarding its sign, is indeed small for niq >400 MeV which 
confirms a consistent trend for the value of niq. Therefore we consider niq =400 MeV to be a 
lower bound for the constituent quark mass in our model. Our results for magnetic moments get 
worse for niq > 500 MeV (see. Table 1). The value of 500 MeV may therefore be considered an 
approximate upper bound for niq. 

Table 2. Contribution of Dirac Form Factor to Nucleon Radii 

comp. Eqs. (19)-(21) 



niq, Ajv in GeV 


F]^{0) in GeV-2 


F5'(0) in GeV-2 


0.315, 2.40 


3.04 


-0.229 


0.320, 2.25 


2.80 


-0.215 


0.330, 2.00 


2.80 


-0.176 


0.340, 1.80 


2.89 


-0.138 


0.350, 1.70 


2.87 


-0.125 


0.360, 1.60 


2.87 


-0.115 


0.370, 1.50 


2.95 


-0.085 


0.380, 1.40 


3.02 


-0.073 


0.390, 1.35 


3.10 


-0.055 


0.400, 1.30 


3.10 


-0.047 


0.410, 1.30 


3.04 


-0.045 


0.420, 1.25 


3.11 


-0.037 


Exp. (in GeV-2) 


2.66 


0.008 



The electric radius of the proton and the magnetic radii of protons and neutrons depend 
fairly weekly on the parameters niq and Ajy. The theoretical values for the neutron electric 
radius squared show the right trend as well, despite of the small size of this quantity. In Table 3 
the best fit of the nucleon static properties for niq =420 MeV and Ajv=1.25 GeV is compared to 
other theoretical approaches such as QCD sum rules [9, 28], the nonlinear chiral model of [4], the 
relativistic invariant quark model in the light-front formalism of [29], the nonrelativistic quark 
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model of [30], the relativistic quark model with chiral symmetry of [13]), the nonrelativistic quark 
model with QCD vacuum induced quark interactions of [31] and the MIT bag model [32]. 



Table 3. Static Properties Compared to Theoretical Approaches 









r^, im 


^ r im 


^ Ml ™ 




Our 


2.79 


-1.86 


0.92 


-0.132 


0.84 


0.84 


loffe [9] 


2.96 


-1.93 










Bahtsky [28] 


2.80 


-1.95 










Meissner [4] 


2.77 


-1.84 


0.97 


-0.25 


0.94 


0.94 


Aznauryan [29] 


2.811 


-1.848 


0.811 


-0.094 


0.825 


0.781 


Isgur [30] 


2.79 


-1.86 










Baric [13] 


2.73 


-1.975 


0.85 








Dorokhov [31] 


2.20 


-1.47 


0.831 









De Grand [32] 


1.90 


-1.27 


0.728 









Exp. 


2.79 


-1.91 


0.86±0.01 


-0.119±0.004 


0.86±0.06 


0.88±0.07 



The shape of the electromagnetic nucleon form factors in the region < < 1 GeV^ for the 
same optimal values of the parameters, mg=420 MeV and Ajv=l-25 GeV, is shown in Figs. 4-7. 
The normalized magnetic form factors of nucleons G^{Q^) / G^{0) are plotted in Fig. 4 and Fig. 5. 

The short dashed line corresponds to the contribution of the triangle vertex function A^^(p,p') 
whereas the medium-long dashed line gives the contribution of the bubble or contact term 

bubbieiP^ P')- This Contribution is roughly constant and not large. It matters however for 
the details of the fitting procedure. The total contribution is marked by the solid line. The 
long-dashed line is the dipole approximation D{Q^) = [1 -|- /0.71GeV^]~^ . It can be taken as 
a representation of the data which it fits quite well. One can see that our results are close to the 
dipole fit and hence to the data. 

The electric form factors of the proton and the neutron are shown in Fig. 6 and 7, respectively. 
The data in Fig. 7 characterize the experimental situation for the neutron electric form factor. We 
remark that the behavior of the nucleon electromagnetic form factors is sensitive to the choice 
of the constituent quark mass, as we have checked. Increasing the constituent quark mass above 
500 MeV and reducing the parameter Ajv below 1.15 GeV improves agreement with the static 
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Fig. 4. Magnetic form factor of proton. Fig. 5. Magnetic form factor of neutron. 
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Fig. 6. Electric form factor of proton. Fig. 7. Electric form factor of neutron. 
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nucleoli properties. However, the form factors at Q^=l GeV^ are then definitely getting too low. 
Thus the value mg=500 MeV is an upper limit in our model also from the point of view of the 
form factors. 

5 Conclusion 

The electromagnetic form factors of nucleons have been calculated within a relativistic three quark 
model with Gaussian shape for the nucleon-quark vertex, and standard (non-confined) quark 
propagators. Gauge invariance of the nonlocal hadron-quark interaction has been implemented 
by the path-independent definition for the derivative of the time-ordering P-exponent. The 
allowed region for the two adjustable parameters, the range parameter Ajv appearing in the 
Gaussian and the constituent quark mass m^, has been obtained by fitting the data for the 
magnetic moments and the electromagnetic radii of the nucleons. The corresponding shapes of 
the electric and magnetic nucleon form factors for mg=420 MeV and Ajv=l-25 GeV are close 
to experiment for < < 1 GeV^. A large body of data is therefore described by only two 
free parameters. It remains to be seen whether a further improvement is possible by introducing 
confined propagators into the model. 

6 Acknowledgements 

We would like to thank P. KroU, R. Rosenfelder, J. Bolz, H. Kamada and H.-Q. Zheng for useful 
discussions. This work was supported in part by the INTAS Grant 94-739 and by the Russian 
Fund of Fundamental Research 94-02-03463-a. 



16 



7 Appendix 



A. Feynman Rules for Nonlocal Electromagnetic Vertices 

We illustrate how to obtain a nonlocal vertex which couples a nucleon, a photon and three 
quarks, see Fig. 2. The analogous nonlocal vertices with mesons were introduced in [23]. Let 
p denote the momentum of the outgoing baryon, k the momentum of the outgoing photon, 
and gi,g2,93 the momenta of incoming quarks, respectively. We start from the gauge invariant 
Lagrangian describing the interaction of nucleons with quarks, Eq.(7). The vertex function 
^^{(liqiqspk) corresponding to Fig. 2 is 

3 

3 (' (' ^(y~] q^Xi—pXi — kxs) 

i{27Ty5(*\p-k-J2q,)Ti:^^^Jq,q2qspk) = j dx,... j dx,e Ri:^^^Jx,...x,){A.l) 

8 = 1 

where the function R'^^^^^^[xi...X5) is given by 



5A^{x5)5N{x4)5q''^{xi)...5q''^-{x3) 



(A.2) 



egNe^^^-^HT^) x (Cr2)F 



xi 



X2 



[tQ] X {t2t) dz^ + TX {Qt2t) dz^ + TX {t2tQ) / dz 



5{X5 - Z) 



We define a typical integral as 



-'typ 



F 



dxi... j dx^e 
which may be written as 

3 

6^3 j dx j d^i j d^2F{^l + ^l)exp ixi^q^ - p) 



(xg - z) 



-'typ 



8 = 1 



(A.3) 



exp 



x-2ii 



^i(-2gi + g2 + gs) + «6V3(g3 - q2) j dz^ exp[ 



-IKZ\ 
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Using the identity 



d^i d^2S{^i)S{^2)F{-dl-dl)m,^2) 



one 



finds 



typ 



3 

8 = 1 



exp[i^i{-2qi + ^2 + gs) + «6^('?3 - q2)] J dz^ e-x^[-ikz\ 

X 

Finally, using the transformations described in Sec. 2 we obtain 



-'typ 



,q2,q3,k) (A. 4) 



8 = 1 



M{qi,q2,q3,k) 



F[A(gi, g2, g3, A:)] - F[A(gi, ^2, g3, 
A(gi, g2, gs, A;) - A(gi, ^2, gs, 0) 



3 3 

A(gi,g2,g3, A;) = + A;^ - ^ g.-g^ + A;(g2 + ^3 - 2gi)]. 

8<i 



8 = 1 



Using the result for the typical integral (A. 4) in Eq.(A.l) we can write the vertex function 
^a^a.aAqiq^qspk) as 



K.a.asiqiq^qspk) = 2iV3eg^e^^^^^^{r^) X (cr2) 



(A.5) 



tQ X (r2f)(g2 + q3- 2gi + k fM{qi, q2, gs, k) 



+ T X {QT2T){q3 + qi- 2q2 + k fM{q2, qs, qi, k) 



+ f X {T2TQ){qi + q2- 2q3 + A;)^M(g3, gi, ^2, k) 
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B. The nucleon mass operator and the jNN vertex 



To demonstrate the calculation of an electromagnetic vertex and nucleon mass operator we 
consider the two generic integrals 



■4 IJ 4 



n 



1 



fB.ll 



F{-I2[kl + k'^ + kEk'E\)F{-l2[kl + k'^ + kEk'^] - + l2kEQ) 



d kjT! f d k rp I , -I — r 

' ^ ' dtJJ 



\ ml + PfE 



F{-l2[kl + k'i + kEk'E\) 



X F'{-I2[kl + k'^ + kEk'^] - t[AQ^ - l2kEQ]) 
with the Euclidean momenta PiE being related to nucleon momenta p^;, by 

Pie = kE + Pe/^, P2E = ^'e + Pe/^, P3E = kE + k'^ - pe/^, 
Pae = kE + p'e - '^Pe/^, Q = Pe - p'e- 



The first integral arises in the calculation of the triangle vertex function or the derivative of the 
nucleon mass operator (in this case pE = p'e) The second integral arises in the calculation of the 
bubble vertex functions. 

00 

Using the Feynman parametrization \ = j c/aexp(— aA) and integrating over kE and k'^ one 



finds for on mass-shell nucleons 



00 00 



/i(-mjv,-mjv, Q ) 







dai...da4 



■ exp 



Ai 



(2 + ai2)(l + 2a34) 



exp 



■((1 + a34)(ai - a2) + (1 + ai)(a34 - "2) + (1 + a2)(a34 - Oi) ) 



3Ai 



exp 



1 00 00 



/2(-m^,-m^,g2) = 12 dt 



000 



dai.-.das 



3QH 



2 + ai2)(2 + 2a3 -t; 
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exp 



^-((1 + a3)(ai - 02)' + (1 + ai)(a3 - 02)' + (1 + a2)(a3 - ai)') 



3A, 



exp 



-12 m 



2x3 



E 

8 = 1 



a,: 



Introducing the notation = ai + we have 

4 33 
Ai = 3 + 2 ^ + ai2a34 + aia2, A2 = 3 + 2 ^ + ^ a^a^. 

i=l i=l i<Cj 

Making a change of variables 

«! = xs, 0,2 = x{l — s), as = yt, a4 = y{l — t) in integral m^, —mif, Q^] 



«! = xs, a2 = x{l — s), as = y, 



in integral /2(— m^, — m^, 



one 



finds 



00 00 1 



/i(-m^, -mjv, 



Q^) = j dx j dy j ds-^exp 



000 



3g; 

' A 



(2 + x)(l + 2y) 



(B.2) 



exp 



^-((1 + y)(l - 2sfx^ + (1 + X5)(y - x{l - sjf + (1 + x{l - s)){y - xsj') 



3A 



exp 



-12(a; + y) m 



2 nij^ 



00 00 1 1 



/2(— m^, — m^, Q^) = 12 J dx J dy J ds J dt-^exp 



0000 



3Q't 



(2 + x)(2 + 2y-t) 



(B.3) 



exp 



^-((1 + y)(l - 2sfx^ + (1 + X5)(y - x{l - sjf + (1 + x(l - s)){y - xsj') 



3A 



exp 



-12(a; + y) m 



2 



" 9 

where A = 3 + 2(a; + y) + a;^s(l — s) + xy. 
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